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Abstract 

The Degasperis-Procesi equation (DP) is one of several equations known 
to model important nonlinear effects such as wave breaking and shock cre- 
ation. It is, however, a special property of the DP equation that these two 
effects can be studied in an explicit way with the help of the multipeakon 
ansatz. In essence this ansatz allows one to model wave breaking as a colli- 
sion of hypothetical particles (peakons and antipeakons) , called henceforth 
collectively multipeakons. It is shown that DP multipeakons have Painleve 
property which implies a universal wave breaking behaviour, that multi- 
peakons can collide only in pairs, and that there are no multiple collisions 
other than, possibly simultaneous, collisions of peakon-antipeakon pairs 
at different locations. Moreover, it is demonstrated that each peakon- 
antipeakon collision results in creation of a shock thus making possible a 
multi shock phenomenon. 

1 Introduction 

The Degasperis-Procesi (DP) equation [9] 

u t - u txx + 4uu x = 3u x ii xx + uu xxx (1.1) 

belongs to a class of one-dimensional wave equations which have attracted con- 
siderable attention over the last decade, following the most studied equation in 
this class, namely, the Camassa-Holm (CH) equation [3] 



Both these equations can be derived from the governing equations under the 
assumption of moderate amplitude [13, 7]. What makes them special is that, 
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on one hand, both are Lax integrable, on the other, both exhibit wave breaking 
phenomenon not captured by linear theory or shallow water, small amplitude 
theory like the Korteweg-deVries equation. The most relevant to this paper 
study of the breakdown of solutions for the CH equation was done by H.P. 
McKean [22, 21]. In these works, it was argued that the breakdown of the CH 
waves is controlled by a kind of caricature of the higher dimensional vorticity, 
namely, m = u — u xx (see [20] ). In particular, it is the initial relative position of 
regions with positive m versus negative m that signals whether the breakdown 
will happen at some later, finite, time. One of the fascinating aspects of both 
the CH and DP equations is the existence of special solutions, pcakons, which 
play the role of basic building blocks of the underlying full theory. Peakons are 
a simple superposition of exponential terms 

n 

u(x,t) = j^ro J -(t)e-l x - x 'WI 

3=1 

for which the function m referred to earlier is m = 2^" =1 nij6 Xj . Were we 
to take the analogy with vorticity at its face value, m for peakons could be 
viewed as a collection of point vortices, situated at Xj(t)s, of strength rrij{t) 
each, initially ordered in some fixed way, say, xi(0) < £2(0) < • • • < x„ (0). The 
case of CH peakons shows that if the strengths {771^(0) : 1 < k < n) arc not of 
the same sign then collisions can occur, meaning that Xj(t c ) = Xj+i(t c ) for some 
j and some time t c . Each collision is accompanied by a blow-up of rrij(t c ) and 
nij + i(t c ) resulting in the derivative u x (t c ) becoming unbounded even though 
u(t c ) remains bounded, in fact continuous. Thus peakons can be used to test 
ideas about wave breaking, the advantage being that the peakon dynamics is 
described by a finite system of ODEs (see Section 2 ). The analysis of the CH 
peakon collisions in this case was done in [1] with the help of explicit formulas. 
In short, the CH peakon problem can be solved by Stieltjcs' formulas involving 
continued fractions [2]. Moreover, the underlying boundary value problem is 
self-adjoint, in fact it is equivalent to an inhomogeneous string which remains 
isospectral under the CH flow. 

The case of the DP equation is superficially similar to the CH case. However, 
deeper analysis shows a remarkable number of new features. For example, the 
associated spectral problem, termed a cubic string in [19], is not self-adjoint 
and this has the immediate consequence that the inverse problem is by far 
more involved. The peakon problem in the case of the positive measure, that 
is when all weights rrijS are positive, was solved explicitly in [19]. However, the 
generalization to the case of a signed measure m is not straightforward since 
the spectral data breaks up into several types depending on the degeneracy of 
the spectrum, as well as on certain coincidental phenomena of anti-resonances 
(eigenvalues Zi, zj pairing according to Zj + Zj = 0). By contrast, the distinction 
between peakons for positive measure m and peakons for the signed measure m 
is less sharp for the CH case where the formulas for peakons can be analytically 
continued from former to latter. This is not so for the DP case. This difficulty 
notwithstanding, in a way analogous to what happens in the CH case, the 
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presence of a collision signals an occurrence of wave breaking; in the DP context 
the connection between wave breaking and peakon collisions was studied earlier 
by H. Lundmark in [17] for the case n = 2 and further by the present authors 
in [23] for n = 3. Important questions of stability and general analytic results 
dealing with DP peakons and the DP wave breaking have been addressed in 
[15, 14, 16, 10]. A considerable amount of work has been also done on adapting 
numerical schemes to deal with the DP equation; we just mention a few: an 
operator splitting method of Feng and Liu [11], or numerical schemes discussed 
by Coclite, Karlsen and Risebro in [6]. 

The DP equation, in contrast to the CH equation, admits shock solutions 
(see [4, 5] for a general, very thorough, discussion). It was H. Lundmark who 
introduced the concept of shockpeakons 

n 

u(x,t) = £{ mj - (t) - Sj (t)sga(x - Xj (t))}e- ]x - Xim , 

3=1 

for which 

n 

m = 2j2{m j 6 Xj +s j S' x .}, (1.2) 

3=1 

and showed that the solution describing a collision of two peakons has a unique 
entropy extension to shockpeakons. He also hypothesized that this might be 
a general phenomenon valid also for n > 2. We prove his conjecture. More 
precisely we prove that the distributional limit of m at the collision point t c 
indeed produces shockpeakon data (1.2) with positive shock strengths Sj thus 
allowing a unique entropy weak extension (see Theorem 5.1 and Corollary 5.2). 

Let us briefly describe our strategy. Instead of analyzing numerous spectral 
types we concentrate on analytic properties of Xj(t),mj(t) as functions of t. 
To this end we analyze the inverse spectral problem for the cubic string with 
the input data of a finite, signed measure. We prove that each Xj(t) must 
be a holomorphic function at t c , while nij(t) is in general only meromorphic 
(Theorem 3.5). Then we perform singularity analysis of the ODEs describing 
peakons (2.6) and prove their Painleve property with the help of Theorems 3.5 
and 3.6 followed by a singularity analysis at the time of collisions described by 
Theorem 4.5. 

The plan of the paper is as follows: we review basic facts about the DP 
equation in Section 2, in Section 3 we discuss the inverse problem for peakons 
of both signs generalizing the uniqueness result known from the pure peakon 
case [19] and use this result to establish analytic properties of positions XjS 
and masses rrijS. In Section 4 we analyze the singular behaviour at the time 
of collisions and establish a universal singularity pattern according to which, in 
the leading term, only the time of the blowup depends on the initial conditions 
not the residue. This fact is proven in Theorem 4.5. We furthermore rule out 
triple collisions in Theorem 4.7, and give an example of a simultaneous collision, 
in different positions, of two peakon-antipeakon pairs; finally in Section 5 we 
prove Theorem 5.1 stating that the distributional limit of colliding peakons is 
indeed a shockpeakon. 
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2 Basic Facts about the DP equation 



The nonlinear equation 

Ut — u xxt + Auu x = 3u x u xx + uu xxx , (2.1) 

often written as 

m t + m x u + 2>mu x = 0, m = u — u xx , (2-2) 

was introduced by Degaspcris and Procesi [9] as an example of a nonlinear 
partial differential equation satisfying asymptotic integrability appearing in the 
family of third order dispersive equations: 

Ut - a 2 u xxt + ju xxx + c u x = (ciu 2 + c 2 u x + c 3 uu xx ) x , 

other examples of integrable equations in this family are the Korteweg-dcVrics 
equation (KdV) and the Camassa-Holm (CH) equation. Formal integrability 
for the DP equation was established by Degasperis, Holm and Hone [8] through 
the construction of a Lax pair and a bi-Hamiltonian structure. In particular, it 
was shown in [8] that the DP equation admits the Lax pair: 

(d x - d xxx )$ = zmf, $t = [z- 1 (l-d*)+u x -ud x }*. (2.3) 

Moreover, one can impose additional boundary conditions provided they do 
not violate the compatibility of these equations. One such a pair of boundary 
conditions was introduced in [19]: 

— > e x , as x -» —oo, $ is bounded as x — > +oo, (2-4) 

where it was also shown that the spectrum of this boundary value problem will 
remain time invariant (isospectral deformation). It suffices for our purposes to 
restrict our attention to the case in which to is a finite discrete (signed) measure. 
Thus for the remainder of the paper we will use the multipeakon ansatz 

n 

u {x,t) = ^ra 1 (t)e- |l - I ' WI (2.5) 

i=l 

where xi(0) < £2(0) < ••• < x„(0) and TOj(O) can have both positive and 
negative values. This ansatz produces m = 2 5Zi=i m i&xi ■ Moreover, with the 
proper interpretation of weak solutions to equation (2.1) we can easily check 
that u is a weak solution to (2.1) provided Xi(t),m,i(t) satisfy the following 
ODEs 

n 

x k (t) = u(x k ) =^ mi (t)e-^- x ^, (2.6a) 

i=l 

71 

rh fc (t) = 2 m k (t) (u x (x k )) = 2 m k (t) £ m^t) sgn(x fc (i) - Xi (t)) e -l**(*)-»*WI , 

i=l 

(2.6b) 
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where (f(x)) = lim e _ ) . + \{f{x + e) + f(x — e)] is the average of / at the point x. 
We will refer to rrijS as masses to emphasize their role in the spectral problem. 
We also need a bit of terminology regarding the phenomenon of breaking. We 
will say that a collision occurred at some time t c if Xi{t c ) = Xi+\(t c ) for some i. 
We can make this concept more geometric by introducing a configuration space 
in which to study pcakon solutions, namely the sector I={x£ R™ \xi < X2 < 
■ ■ ■ < x n }. Then a collision corresponds to the solution Xi hitting the boundary 
of X. 

A very useful property of equations (2.6) is the existence of n constants of 
motion. This follows readily from Theorem 2.10 in [19]. 

Lemma 2.1. M p (1 < p < n), given by: 

m p= e (ii m ) fn(i-e^-^) 2 ] 

( Je £- fcl ) V*ei / \j=i J 

are n constants of motion of the system of equations (2.6) ; where ('■ 1 ' n '\ is the 
set of all p-element subsets I = {ii < • • • < i p } of {1, . . . , n}. 



3 Inverse Problem for multipeakons 

The boundary value value problem (2.3) and (2.4) can be transformed to a finite 
interval boundary value problem, the cubic string problem. Indeed, following 
[19], the change of variables (Liouville transformation) 



y = tanhf, *(x) = ^4 (3- 1 ) 

z i- — y 



maps the DP spectral problem into the cubic string problem: 

-<t>yyy{y) = zg{y)<j){y), -1< y < 1, (3.2a) 
0(-l)=^,(-l) = 0(l)=O, (3.2b) 

where g is the transformation of the measure m induced by the Liouville trans- 
formation (3.1). Furthermore, as one can explicitly check, g is also a finite 
signed measure and its support does not include the endpoints if the original 
measure to is a finite signed measure. More concretely, in this paper, 

II: 

.9 = E 5l<5 ^' - l <Vi <V2 < ■■■ <Vn <l, (3.3) 

1=1 

with weights gi € R. The inverse problem is studied with the help of two Weyl 
functions. 
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Definition 3.1. Let <fi(y;z) denote the solution to the initial value problem 
(3.2a) with initial conditions <f>(— l;z) = <j> y (—l;z) = 0, (f> yy (—l;z) = 1. The 
Wcyl functions are ratios: 

W(z) = ^, Z ^-^TT- 

<j>{l;z) nl;z) 

These two functions encode spectral information needed to solve the inverse 
problem. It is easy to verify that in the case of (3.3) both W(z) and Z{z) are 
rational functions which makes inversion algebraic. However, in contrast to the 
pure peakon case gi > 0, the spectrum of the boundary value problem (3.2a) is 
in general complex and not necessarily simple. This makes the inversion more 
challenging. Regardless of the complexity of the spectrum though the Weyl 
functions undergo a simple evolution under the DP flow. Indeed, using the 
second member of the Lax pair given by (2.3) one can find the time evolution 
of W{z) and Z[z). To wit, using results from Theorem 2.3 in [23] we obtain the 
following characterization of the time evolution of W(z) and Z(z). 



Theorem 3.2. Let 



be the partial fraction decomposition of w ^ z > } where dj denotes the algebraic 
degeneracy of the j-th eigenvalue. Then the DP time evolution implies: 

(V 

bf =pf\t)e$, 
where Pj (t) is a polynomial in t of degree dj — k. 

(2) m+^eLi^^E,^- 

(3) 

W - 1 

W= + M+, Z = (W - 1)M+ + W 

z 

We immediately have: 

Corollary 3.3. Under the DP flow the Weyl functions W, Z are entire functions 
of time. 

The uniqueness result below plays a major role in the solution to the inverse 
problem. 

Theorem 3.4. Suppose $ : g — > {W(z), Z(z)} is the map that associates to the 
cubic string problem (3.2) with a finite signed measure g, the Weyl functions 
W(z),Z(z). Then $ is infective. 



G 



Proof. The proof relies on remarks made in [18]. We will construct a recursive 
scheme to solve the inverse spectral problem; given W and Z obtained from the 
map $ we will reconstruct the finite, signed measure g whose Weyl functions 
are W and Z. More precisely, we will show that the yj's and g/s in equation 
(3.3) are uniquely determined from W(z), Z(z). First we recall that W and Z 
are constructed from solutions to the initial value problem (sec Definition 3.1) 

-^yyyiv) = zg{y)4>{y), —l<y<l, (3.4a) 
0(-l)=^,(-l) = O, <M-1) = 1 (3.4b) 

Masses gj are situated at yj, 1 < j < n and for convenience let us set yo = 
0, j/n+i = 1 and denote by lj = yj+i — yj the length of the interval (yj,yj+i). 
Then on each interval (yj,yj+i) the solution to (3.4) takes the form 

<fi(y) = 0(%+i) + Mvj+i) (y - Vj+i) + ^wivj+i-) (y - y j+ i) 2 /2, o<j<n 

and the condition of crossing yj+i is: continuity of <f> and <j> y and the jump 
condition 4> yy {yj+i+) — 4>yy{yj+i — ) = ~ z 9j+i4>{yj+i)- We establish, for example 
by an easy induction, 

(0( %+1 ),^(j/ J+1 ),^(y J+1 -)) = {- Z y n (l*/2,l j: l) +0(z j ~ 1 ), 

k=l 

(3.5) 

valid for < j < n, with the convention that for j = the product equals 1 
and there is no remainder. Likewise, 

j ad 2 

ttfa)MVi),<ihn l (Vi+)) = LI (0,0,l) + O(^- 1 ), (3.6) 

fc=l z 

valid for 1 < j < n. 

For < j < n we define (w 2 j,z 2 j) = (^-, ^f-)\y= Vj+1 - and (w 2 j-i, z 2 j-i) = 
(^"' ljr~)\v=yj+- These quantities are essentially the left hand and the right 
hand analogs of Weyl functions introduced in Definition 3.1 and correspond to 
shorter strings terminating at yj+i with no mass at the endpoint, or terminating 
at yj but with the mass gj at the end. Equation (3.4) implies that the sequence 
(w 2 j 1 Z2j, v)2 j—i, Z2j-\) satisfies the recurrence relations 



W2 3 7 1 7 W2 3 , "3 I-} 1\ 

w 2j -i = — J --lj, z 2j - 1 = — -l j —± + f, (3.7) 

Z 2 j z 2] z 2] I 

w 2 j-2 = , z 2j -2 = V zg 3 \ (3.8) 

Z 2 j-1 Z 2 j-1 

the iteration starts at W2n = W(z),z 2n = Z(z) and terminates at W—i,Z—\. 
Moreover, based on equations (3.5) and (3.6), we easily establish 
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2 12 1 1 

W 2 j = — + 0(-), Z 23 = j2+0(-), W2j-1 = Z2j-1 = O(-), as z oo, 

Lj Z Lj Z Z 

(3.9) 

which implies that the quantities {lj,gj} are determined in each step from the 
large z asymptotics of terms known from the previous step. Indeed, if we denote 
by a' 1 ™* 1 the coefficient of z~ m in the expansion of a holomorphic function a(z) 
at z — oo we obtain the recovery formulas 

2 1 
l i = — foT' 9j = — 7u- ■ (3-10) 

Thus we proved that given a pair of Weyl functions W(z), Z(z) obtained from a 
cubic string problem (3.4) with a finite, signed measure g, there exists a unique 
solution to the recurrence relations (3.7) subject to (3.9) and thus a unique 
cubic string corresponding to W(z), Z(z). □ 

We are now ready to state the central theorem of this section 

Theorem 3.5. Let {xj(t),mj(t)}, j = 1, ...,n be the positions and masses 
of the peakon ansatz (2.5) corresponding to an arbitrary signed measure m = 
2 y~]j— i rnjS Xj , satisfying peakon equations (2.6) on the time interval (0,t c ) and 
suppose that a collision occurs at t c . Then the positions x±(t) . . . , x n (t) are ana- 
lytic functions at t c , while the masses m\{t) . . . m n (t) are given by meromorphic 
functions at t c . 

Proof. Given the initial conditions {xi(0) < £2(0) < • • • < x n (0)} and 
{mi(0), 7712(0), . . . ,m ra (0)} we set up the string problem (3.2a) after mapping 
m(0) to ,g(0). This produces the Weyl functions W(0), Z(0), which under the 
peakon flow evolve as entire functions of time in view of Corollary 3.3. We then 
set up the recursive scheme (3.7) with W(t),Z(t) as inputs. At each stage of 
recursion only rational operations are involved and since the recursion is finite 
the formulas (3.10) result in functions meromorphic in t. Thus all gj,yj are 
meromorphic in t. For t < t c all distances lj > and at t c some vanishes but 
all lj remain finite, because this is a finite string. Hence lj(t) is regular at t c 
hence analytic there. For a signed measure g there are no bounds restrictions 
on individual gj so in general gj remains meromorphic at t c . Mapping back to 
the real axis is afforded by y = tanh f ; hence positions of individual masses are 
given by Xj = In y3 ~^ > L . The only singular points of this map are for yj = ±1 
which means the end of the string or, after mapping the problem back to the 
real axis, ±00. However, based on results in [23], none of the masses can escape 

p.. y ■ -1-1 

to ±00 in finite time. So J i—1 is in the domain of analyticity of In and hence the 
XjS are analytic at t c . The relation between the measures m and g appearing in 

equations (2.3) and (3.2a) is given by rrij = — ^ — gj which implies the claim 
since gj is meromorphic and yj analytic. □ 



The above theorem establishes that the only singular points of solutions to 
the peakon ODE system (2.6a) and (2.6b) are poles. Since the inverse problem 
argument is valid for a fixed ordering x\ < xi < ■ ■ ■ < x n of masses, the 
analytic continuation of masses and positions into the complex domain in t will 
satisfy equations (2.6a) and (2.6b) in which sgn{xu — xi), e~\( Xk ~ Xi >> are replaced 
with sgn(fc — i), e~ s e n ( k - l ){ x k-xi) respectively, to be consistent with the original 
ordering. It is for these equations that we note the absence of movable critical 
points also known as Painleve property [12]. To facilitate the statement of the 
last theorem of this section we set Xi = e Xi , 1 < i < n and rewrite the system 
(2.6a) and (2.6b) in new variables {nii,Xi}. 

Theorem 3.6 (Painleve property). The system of differential equations 

n ( X \ sgn ( fe ~*' ™ / y. \ s s n ( fe — 

Xk = X k mi f j i ^ = 2111^™, sgn(fc -i) ij^J 

has the Painleve property. 

Proof. First we observe (using the variables of the proof of Theorem 3.5) that 
X* = hence X^s arc meromorphic in t because so are j/jS. The formulas 

for XiS and to^s obtained from the inverse problem are meromorphic in t in 
the complex plane C and depend on 2n constants (spectral data consisting, in 
the generic case, of n positions of poles and n residues of the Weyl function 
W), which for the cubic string problem, in view of the ordering condition, 
are confined to an open set in C 2n by continuity of the inverse spectral map. 
Relaxing that condition results in a solution depending on 2n arbitrary constants 
which comprises a general solution which is meromorphic in t in the whole 
complex plane C. □ 

In the remainder of the paper we will concentrate on the specific singularity 
structure at the time of collisions of peakons. 



4 Blow-up behaviour 

We now proceed to establish several theorems on peakon collisions for DP equa- 
tion. To begin with we recall the definition of a peakon collision briefly discussed 
in the introduction. We call t c the collision time if there exists some i such that 

lim Xi(t) = lim Xi+i(t), (4-1) 

where Xi(t)s are the position functions in the ansatz (2.5). Equivalently, we say 
that the z-th peakon collides with the (i + l)-th peakon at the time t c . If there 
exist more than two position functions being identical at t c then we will say 
that a multiple collision happens at t c . 

In this section, wc describe the behaviour of the peakon dynamical system 
(2.6) in the neighbourhood of a collision time t c . 



9 



To this end we need to study a special skew-symmetric nx n real matrix A n 
given by 

A n = [sgn(i - j)aij] (4.2) 
whose entries satisfy = aji and 

fly = o-ao-ij 7^ 0, for all < i < I < j < n. (4-3) 

The following propositions hold for such a matrix. 



Lemma 4.1. There exists a matrix P with detP = 1 such that P T A n P = B n , 
where 



-<2i2 

a\2 



-a 34 
a 34 



— a n _i, 



i/ n is ewn, 



or 



S„ = 





&23 



-a 2 3 






045 



-0,45 




, i/ n is odd. 



-In-Ill 



Proof. The conclusion is trivial for n = 1,2. We assume the conclusion to hold 
for n — 2; to show that it holds for n we divide A„ into four block submatriccs 



A, 



-fln-ln 
dn-ln 

direct computation shows that 



where C 



A„_ 2 — -B 
Let us set Pi 



In-2 

-C- l B T h 



P(A n P x 



A n _ 2 - BC- X B T 
C 



(4.4) 
, then a 



In view of condition (4.3), B can be written as (a, a„_i „a), where a = {a\ „_i, 0,2 n-i; ■ • ■ j ftn-2 n-i) T - 
It is now elementary to verify that BG~ X B T = 0. By the induction hypothesis 



10 



there exists a matrix P 2 with detP2 = 1 such that Pj A n _ 2 P 2 = B n -2, hence if 



we set 



the conclusion follows. 



P = Pi 



P 2 
I 2 



□ 



Corollary 4.2. If n = 2k then dct A 2k = Y[ a 2i-i 2 i > °- Ifn = 2k + 1 then 
the rank of A 2 k+i is 2k. 



Lemma 4.3. Let E = (1, 1, . . . , 1), n = 2k + 1 and all entries satisfy < Ojj < 

E 
A 2 u+i 



1 . Then the rank of the matrix 



is 2k + 1, 



Proof. Let n be any odd number. It suffices to show that the determinant of 
the matrix 



A n — 



1 


1 


1 




1 


012 





-a 2 3 


-a 24 


— a 2n 


013 


a 2 3 





-a 34 


~03„ 






034 







• '■ '■ — 0"a-1n-\ —O-n-ln 

■ a>n—2n— 1 — a n —i n 

Oln 02 n 13 n •'• •*• On-ln 

is positive. For n = 3 direct computation shows that 

det j4 3 = ai 3 (l - a 23 ) + a| 3 > 0. 

We assume now that the conclusion holds for n — 2. We will show that it also 
holds for n. First, we divide A n into four submatrices by 



where 

C = 





On — 1 1 



"O n _i, 





B = 



B T 



-B 
C 



Oln-l Oln 

b a„_i„b 



B = 



-1 -1 
b a„_i„b 



and b = (a 2n ^\ 1 a 3 „_i, . . . ,a„_2n-i) T - Since C is invcrtible we can factor A Tl 
into the product of upper and lower block triangular matrices as follows: 



A n - 2 -B 
B T C 



A n _ 2 + BC- x B T -BC- 1 
h 



I n - 2 
B T C 
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Hence deti„ = det(i„_ 2 + BC- l B T ) detC = a 2 n _ ln det(I„_ 2 + BC- l B T ). 

Direct computation shows that all the entries of BC~ 1 B T vanish except the 
first rowwhich equals (a~_ ln — l)(ai n _i, a2 n -ii • • • ; o,n-2n-i), therefore 



det(i„_ 2 + BC- l B T ) 

= dct 1„_ 2 + (K-ln - 



ai2 



ai3 



023 



Oln-2 d2n-2 



-a 2 3 




-034 



On-2n-l 
— fl2n-2 
— 03„_2 



fln-3n-2 



deti„_ 2 + (a n _ ln - 1) 



Oi2 — 023 •■• ••■ — &2n-2 

ai3 a 2 3 — a 3 4 •■■ — a 3 „_ 2 

Oln-2 fl2n-2 ,-- •*• 0, n -3n-2 

flln-1 fl2n-l ••• ••' fln-3n-l 1n-2n-l 



where we used that n is odd. Finally, in view of equation (4.3), we can replace the 
matrix in the second determinant by an upper triangular matrix by performing 
appropriate column additions, obtaining 



det An 



l n — l,n 



deti„_ 2 + K_i„ - 1) 



012 



a 2 3 







det A„_ 2 



«n-2n-l 

l)ai 2 a 2 3 • • • a„_2n-i > 0. 



□ 



By using the lemmas above, we can obtain the property of rrii(t) at the time 
of blow-up. 

Theorem 4.4. If mi blows up at some to then rrii has a pole of order 1 at to. 

Proof. Since TOj's are meromorphic in t we can assume that the leading term 

Ci 7^ 0. If the conclusion does 



in the Laurent series of rrij around to is 
not hold then 



C, 



(t-t ) a < ' 

a = niaxja;} > 2. 



Set S = {ij : = a} = . . . , £&} where ix < i<i < ■ • ■ < and k is at 
least 2 by virtue of Lemma 2.1 with p = 1. Comparing the leading term of both 
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sides of (2.6b) with ij 6 S, one can see the leading term on the left hand side 
is (f-t^o+i while the leading term on the right hand side is 



2C, 



s s n (* 



-Xi,\Q. 



1=1 



Since 2a > a + 1, the coefficient of (t — to)~ 2a must be zero, which leads to a 
homogeneous linear equations A^C = 0, where Ak = (sgn(j — l)aji) is a fc x fc 
skew-symmetric matrix with aji = e x ^ ~ x '' (1 < j < I) and C — {Ci 1 , . . . , Ci k ) T . 
Additionally one can also find 



Cix + Ci 2 + ' ' ' + Ci k — 

by comparing the leading term in Mi. It is clear that Ak satisfies (4.3) and the 
condition in Lemma 4.3. Hence C must be zero according to Corollary 4.2 and 
Lemma 4.3, which leads to a contradiction. □ 

In the proof above, we only use that the rrii(t)s are meromorphic. However, 
we can get stronger conclusions if we also take into account that the Xi(t)s are 
holomorphic. 

Theorem 4.5. // rrij 1 , . . . , rrij k (1 < ji < ■ • ■ < jk < n) blow up at t c and all 
other rrii remain bounded, then the following conclusions hold. 

(1) k must be even. 

(2) t c must be a collision time. Moreover for all odd I such that 1 < I < k, 
the peakon with label ji must collide with the peakon with label ji+i . 

(3) The leading term of m Js (t) in the Laurent series around t c must have 



the form 



(~i) 3 

2(t-t c ) 



for all 1 < s < k. 



Proof. Assume that rrij 1 , nij 2 , . . . , nij k blow up at t c . Since Mi = mi + mi + 
• • • + m n is conserved, k is at least 2. Moreover, by Theorem 4.4 the leading 
term in each 's Laurent series has the form -r^r- Hence, by equations (2.6b), 
the coefficients C = (Ci, . . . ,Ck) T satisfy the linear equations 



A k C = -hi, 



A) 1 



(4.5) 



where the matrix A k = (A\^)i<i tm < k = (sgn(^-m)e Xl ^)i<i, m <k satisfies 

(4.2) and (4.3). Likewise, comparing the leading terms of both sides of (2.6a) 
with the subscript j s (1 < s < k), one finds 



B {k) C = 0, 



(4.6) 



where the matrix = )i<i,m<fc = (e ^ Xjm Xj > \)i<i, m <k- Now we prove 
that the theorem holds for k = 2. In this case, (4.5) and (4.6) reduce to 






-a-12 




' Ci ' 


1 


' 1 " 




1 a 12 




' Ci ' 


ai2 







c 2 


_ _ 2 


1 




ai2 1 




c 2 



0, 
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where ayz = e 1 " x 
tions exists iff a\i 



Direct computation shows that the solution of the equa- 
1, and the solution is C2 = —C\ = \. Since ari = 1 



is equivalent to Xj 1 (t c ) — Xj 2 (t c ), we conclude that the peakon with label j\ 
collides at t c with the peakon with label 32- Suppose now the conclusions are 
valid for k — 2. We will show that they hold for k as well. Let us use the same 
block decomposition as in equation (4.4), obtaining 



-4, 



A k -2 
T 

a 1 

ak-ika 1 







-a/b-i fe 




B (k-2) 
T 

a 1 



a a fc _i fe a 



where a = (aifc_i, a 2 k-i, ■ ■ - ,afe-2fe-i ) T . Let us now combine the last two rows 
of (4.5) and (4.6), writing them collectively as 





a T 





— ftfc-lfc 




1 - 

1 




T 

-lk& 
'I' 

a 1 


fe 

1 




Ofc-1 k 


c = 


2 






T 

-lk& 


flfe-1 fe 


1 








The latter expression can subsequently be easily reduced to 









Ci 






r t 
a 1 





— fe 








-1 





flfc-lfc 


a fe-l fe 




Cfe-2 


1 


-(1 - afc-ifc) 





1 


2afc-i fe 




_ 2 


1 








1 




Cfe-i 
Cfe 




1 



which implies the condition for the existence of the solution to be a\_ lk 
hence afc-i, fe = 1. The latter condition indicates the collision of m,j k with rr 
Futhcrmore, the solution for the last two components of C is then Ck-i ~- 
and Ck — 5, which proves the sign statement for the last two compor 
Substituting ak~ik = l,Cfe_i = — \,Ck = \ into (4.5) and (4.6) and den 
the first k — 2 components of C by C we obtain the following equations: 



1 



^fc- 2 C = — (1, . . . , 1) T , B^-^C = 0, a T C = 0. 



The first two equations hold by the induction hypothesis. To show that the 
third equation holds automatically if the induction hypothesis is satisfied, we 
observe that as the result of collisions (jith mass collides with j2th mass etc.) 
a T = (aifc_i,ai fc _i,a3fc-i,a3fc-i,- • • ,a k -3k-i,o,k-3k-i), hence, indeed, the 
last equation follows from the induction hypothesis. 

□ 
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The following amplification of item 2 in the above theorem is automatic. 

Corollary 4.6. Suppose rrij 1 rrij k (1 < j\ < ■ ■ ■ < j k < n) blow up at t c and 
all other rrii remain bounded. Then for all odd I such that 1 < I < k, the peakon 
with label ji must collide at t c with the peakon with label + 1 ( its neighbour to 
the right). 

So far we have established that when the masses become unbounded the 
collisions must occur. The converse turns out to be valid as well. 

Theorem 4.7. For all initial conditions for which M n ^ the following prop- 
erties are valid: 

(1) If the i-th peakon collides with another peakon /peakons at t c , m^t) must 
blow up at t c . 

(2) For all i, rrii(t) cannot change its sign. In particular, m,i(t) ^ fort < t c . 

(3) There are no multiple collisions. 

(4) The distance between colliding peakons has a simple zero at t c . 

Proof. Without loss of generality we can suppose that the labelling is chosen 
so that xi{t c ) = ■■■ = x kl (t c ) < x kl+1 (t c ) = x kl+2 {t c ) = ••■ = x k2 (t c ) < • • • < 
Xki(tc) for all colliding peakons at distinct positions x kl (t c ) < Xk 2 (t c ) < ■ ■ ■ < 
x kl (t c ). Let us now denote the set indexing all colliding peakons by /. Since 
M n = mi?ri2(l — e Xl ~ X2 ) 2 ■ ■ ■ m„_im„(l — e 1 "- 1 " 1 ") 2 is conserved and nonzero, 
it is clear that some of the masses must become unbounded. Let us denote the 
set of labels of those masses which blow up at t c by J. By Theorem 4.4 any 
such mass corresponds to a colliding peakon; thus J G I. Moreover, any such a 
mass has a simple pole at t c . On the other hand, for each colliding peakons with 
adjacent indices j and j + 1, t c is a zero of (1 — e Xj ~ Xj+1 ) 2 of order bounded from 
below by 2. Thus the order of the zero of all such exponential factors appearing 

in M n is bounded from below by 2([fci - 1] + [ft 2 - 1] H [h - 1]) = 2(|J| - I), 

where |/| denotes the cardinality of /. Hence, since all unbounded masses have 
poles of order 1, \J\ > 2(|J| — i) to ensure that M n ^ 0. The maximum of I occurs 
when the masses collide in pairs, hence I < ^ and thus |J| > 2(\I\ - J|l) = |/|. 
This proves that J = I since J C I and thus (1) is proven. To prove (3) we 
return to the inequality above which now reads |/| > 2(\I\ — I), implying I > 
Since the right hand side is the maximum of I, I = ^ follows, which in turn 
implies that all collisions occur in pairs, hence absence of multiple collisions. To 
prove (4) we note that for M n to remain bounded the order of the zero of all 
exponential factors has to be exactly |/| = 2^- hence each factor (1 — e Xj ~ x i+ 1 ) 2 
has zero of order exactly equal 2. 

This concluded the proof of (1), (3) and (4). In order to prove (2) we 
suppose that for some i, mi(t) changes its sign, then there exists some to for 
which mj(to) = while all m,jS remain bounded since to < t c . Hence 

M„ = lim M n = lim mim 2 • ••m„(l - e Xl ~ X2 ) 2 ■ ■ • (1 - e *n-i-* n \2 = . 

t->i Q t-t-to 
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This contradicts M n ^0. □ 

Remark 4.8. It is now not difficult to verify that the constants of motion 
Mi, . . . , M n can be extended up until the collision time t c by using Lemma 2.1, 
followed by Theorems 4.5, 4.7. 

Corollary 4.9. // mj and rrij + i collide at t c > then mj > and mj+i < 
before the collision. 

Proof. Since collisions only occur in pairs, the leading terms in nij and rrij+i's 
Laurent series must be T 2 (t-t ) respectively. This implies 

lim nij = +oo, lim m J+ i = — oo. 

The conclusion holds since in view of Theorem 4.7 mj and mj+i cannot change 
their signs. □ 

The following proposition shows that the simultaneous collisions (several 
peakon-antipeakon pairs collide at distinct locations at the common time t c ) 
can happen. We indicate below how certain symmetric initial conditions will 
lead to simultaneous collisions. To this end we consider equations (2.6) for n = 4 
and a special choice of initial conditions. 

Lemma 4.10. If the initial conditions satisfy 

toi(O) = -m 4 (0) > 0, m 2 (0) = -m 3 (0) < 0, (4.7a) 
si(0) - x 2 (0) = ar 8 (0) - aa(0) < 0, (4.7b) 

then mx(t) = —n%i(t), m^it) = —m^(t), x±(t) — X2(t) = x$(t) — x^(t) will hold 
for all < t < t c . 

Proof. Consider the following ODEs 

xi = mi + m 2 e 1 1 — m 2 e 1 s — mie , 

x 2 = mie Xl ~ X2 + m 2 - m 2 e X2 ~ X3 - mie 11 " 13 , 

x 3 = mie Xl ~ X3 + m 2 e X2 ~ X3 - m 2 - mie 11 " 12 , 

mi = -2mi(-m 2 e Xl - X2 +m 2 e Xl - X3 +mie 2xi -' X2 - X3 ), 

m 2 = -2m 2 (mie Xl ~ X2 + m 2 e 12 ^ 3 + mie Xl ~ X3 ), 

then direct computation shows that {xi, x 2 , X3, X2+X3 — xi, mi, m 2 , —m 2 , — mi} 
satisfy the system of ODEs (2.6) for n = 4. □ 

The following is then immediate (see figure 1). 

Corollary 4.11. If the initial conditions (4-7) hold and the peakon-antipeakon 
pair (mi,m 2 ) collides at t c then so does (7713,7714) and vice-versa. 
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0.1 0.2 0.3 0.4 O 0.1 0.2 0.3 0.4 

Figure 1: Two symmetric peakon-antipeakon pairs with masses mi(0) = 10 = 
—7714(0), 7712(0) = —1 = —7773(0) undergo a simultaneous collision. 

5 Collisions and shocks 

In this section we investigate the behaviour of 777 and u at the time of collision(s). 
We start with m and observe that since the collision of peakons occurs in pairs 
it is sufficient to study a fixed colliding pair mj,m,j+i. 

Theorem 5.1. If mj collides with rrij + i at time t c > and the position x c , 
then 

lim (nij(t)8(x — Xj(t)) + nij + i(t)S(x — Xj + i(t))) 

t->tc 

= ( lim (mj + 777 3+ i) ) 6(x - x c ) + - (xj(t c ) - x j+1 (t c j) S'(x - x c ) 

= ( lim (irij + 777j + i) J 5(x — x c ) + — ( lim (u(xj(t),t) — u(xj+i(t),t) J S'(x — x c ) 



Proof. For an arbitrary tp(x) € 
(mj (t)5(a;— (t))+m i+ i (i)<J(x-Xj + i (i)), y>(x)) = mj (t)<p(xj (t))+m J+1 (t)ip(xj + i (t)). 
Using Corollary 4.9 we can write 

+ C* + 0(t - t c ), TOj+1 = ~~~r + C* + 0(t - t c ) 



2(t-t c ) u v L; ' j+1 2(i-t c ) 
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around t c . Hence 

lim (rrij{t)5(x — xj(t)) + rrij + i(t)5(x — Xj+i(t)), tp(x)) 



--(C + Co)f(x c ) - lim 



t->tc 2(t-t c ) 
{ lirn (rrij + m J+ i)J <p(x c ) - ^ ( lirn (±j - x j+ i)J <p'(x c ) 

lim (rrij +m j+1 ) J </j(x c ) - - ( lim (u(xj(i),t) - u(x j+ i(t), t)) J y'(x c ), 



where in the last step we used equation (2.6a). The conclusion now follows from 
the definition of (5 and <$' . □ 

Since m = u — u xx we have the immediate corollary. 

Corollary 5.2. Suppose m(0) = 2 ^fe=i Tnk(0)S(x — x k (0)) is a multipeakon 
at t = for which M n ^ and such that at t c one, or several of its peakon- 
antipeakon pairs collide. For any colliding pair k, k+1 let us denote lim^j- u(xk(t)) = 
ui(x k (t c )), lim t ^ t - u(xk+i(t) = u r (xk(t c )) respectively. Then 

n 

hm m(t) =2 J"m k (t c )S(x-x k (t c ))+2 V s k (t c )S'(x-x k (t c )) inV(R). 
t—^t~ — — 

fc— 1 k:pairs 

x k collide 

The shock strengths are given by 

Ul(x k {t c )) ~ U r {x k {t c )) 

Sk{t c ) = g ' 

and they satisfy the (strict) entropy condition s k {t c ) > 0. 

Proof. It suffices to prove the claim if there is only one colliding pair; the general 
case follows easily since masses collide pairwise. For t < t c the measure evolves 
as m(t) = 2 J2 k =i m k{t)S( x ~ x k (t)) where x k (t), m k (t) satisfy equations (2.6a), 
(2.6b) respectively. Suppose now that the pair j,j + 1 collides at the point 
x c . Then by Theorem 5.1 lim t _^ c m(t) = 2 J2k^ 3 m k (t c )5(x - x k (t c )) + 
2 lim t ^ t - (rrij +irij + i(t)5(x — x c ) +2^(xj(t c ) — Xj+i(t c ))6' (x — x c ). To prove that 
s j(tc) > we write Sj(t c ) = h^ m t^t-(xj(t)—Xj+i(t)) = ^lim t _^ t -(u(xj(i),t) — 
u(xj+\(t),t)) and observe 



xj(t c ) - x j+ i(t c ) = lim 



X j+ i(t) - Xj(t) 



tr-t 



which implies the entropy condition Sj(t c ) > in view of the ordering assump- 
tion Xj(t) < Xj + i(t). The strict inequality follows from item (4) in Theorem 
4.7. □ 
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The following amplification of the previous theorem brings the issues of the 
wave breakdown and a shock creation sharply into focus. To put our result 
into the proper perspective we first review the wcll-posedness result for L 1 (M) n 
BV(M) proven by Coclite and Karlsen ([4], Section 3). We present only the core 
result pertinent to our paper. 

Theorem 5.3 (Coclitc-Karlscn). Let u G L 1 (R) DBV{M). Then there exists 
a unique entropy weak solution to the Cauchy problem u\ t= o = uq for the DP 
equation (2.1). 

It is then proven in [17] that the shockpcakon ansatz 

n 

u(x,t) = YsMt) ~ Sj(t) S gn(x - ^(i))}e-l*-^«l, 

3=1 

is an entropy weak solution provided the shock strengths Sj > 0. This sets the 
stage for the next theorem. 

Theorem 5.4. Assume that a multipeakon solution u(x, t) exists on W 1 x [0, t c ), 
then u(;t) £ L 1 (R) n BV(R) for all < t < t c and u(-,t) converges in L 1 to 
the shockpeakon 

n n 

u(x, t c ) = Y, rn t {t c )e-\ x - x ^ + ^ CiXi(t c ) S gn(x - x t (t c ))e^ x ~ x ^ , 

i=l i=l 

u(;t c ) G L^M) n BV(R), where m 4 (<) 's Laurent expansion around t c is written 
as 

mi{t) = — l — +53oi(i-i c ) i = -j-^y+nnit), 

c 1=0 c 

with the proviso that Ci = if the ith mass is not involved in a collision and 
Ci — — | for a colliding peakon, Ci = h for a colliding antipeakon, respectively. 

Proof. We start with the case n = 2. Then u(x,t) = mi(f)e"' I_I1 ' f " + 
^We - ' 1-12 ''" and Xi(t c ) = X2(t c ) — x c . According to Theorem 4.5, we 
have 

™i(*) = -^y +™i(*)> "»>(*) =2(4^+^*)' 
where fh\{t) , rri2{t) are analytic around t c . It is clear that 

mi^e-^-^W 1 + fh 2 (t)e-\ x - X2 ^ G L 1 (R) n BV(R) for all < t < t c . 
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By the mean value theorem we find that 

dcf 1 



v(x, t) 



\x-x 2 (t)\ _ e -\x-xi{t)\\ 



( 1 

2 

1 

2 



2(i - t c ) 
ii(s)e a; - :cl W - x 2 {s)e x - X2[s) 

1 



s =t+e 1 (t c -t) 



2(t - t c ) 



s=t+e 3 (t c -t) 
(i 1 (s)e Xl(s )- :r + x 2 {s)e x - X2[s) ) 



x < Xl(t) < x 2 {t), 
x\(t) < x 2 (t) < x, 
Xi(t) < x < x 2 (t), 



s=t+e 3 (t c -t) 



where < 8j < 1, j = 1, 2, 3. Hence we have the pointwisc limit 



lim v(x, t) 
Let us define 
v{x,t c ) = 



sgn(x - x c )(-\xi{t c ) + \x 2 {t c ))e l x x ^, for x ^ z c , 



Ai(t c )+a;2(tc) 
2 



, for a; = x c . 



sgn(x - x c )(—^xi(t c ) + ±x 2 (t c ))e |x Xc| , for x ^ x c 
0, for x = x c . 



and consider the integral 

+oo 

i) — v{x, t c )\dx 




\v(x,t) — v(x,t c )\dx. 



Then the first and the last term of the right hand side converge to zero as t — > t' c 
due to Lebesgue's dominated convergence theorem. 
Observe that the second term satisfies 

x 2 {t) x 2 (t) x 2 (t) 

\v(x, t) — v(x, t c )\dx < / \v(x,t)\dx + / \v(x, t c )\dx 



Xl(t) 

x 2 {t) 



x 2 {t) 



< 



\v(x, t)\dx 



xi{t) 
x 2 (t) 



xi(t) 



Xl(t) 

,x-x 2 (t) _ e x c -x^ 

2(t c - 1) 



Xl(t) 



x 2 {t) 



x 2 (t) 



d^+o / \xi(s)\e xi ^~ x dx + - / \x 2 (s)\e x - X2{s) dx 



xi(t) 



xi(t) 



xi(t) 



\v{x, t)\dx + X f { ~jf \e«~ x ° - e x °~y\ + \ 



2(0 



■An 



\x 1 (s)\e Xl ^- x dx + - j \x 2 {s)\e x - X2 ^dx 



xi(t) 



Xi(t) 



where s e (t,t c ) and y 6 (x 1 (t),x 2 (t)). Since \xi(s)\e Xl ^~ x and \x 2 (s)\e x ~ X2 ^ 
are bounded, and 

X2(t)-n(t)-»0, ^§TZ^T -»gi(*c)-A 2 (t c ), lew-*" - e*«-»| -> 
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as t — > t~ , we have that v(x, t) converges to v(x,t c ) in the sense of L 1 , which 
shows that the conclusion holds for n = 2. 

In general, since collisions can only occur in pairs, we can assume that 
m,j 1 (t), m J1+ i(t), m j2 (t), m j2+ i(i), . . . , rrij k (t), m, Jk+ i{t) blow up at t c and all the 
other mi's remain bounded. It is clear that mi(t)e~\ x ~ Xi ( t >\ lies in L X (R) n 
BV(M) and converges to rrii(t c )e~^ x ~ Xi ^ tc ^ in L 1 if rrii(t) remains bounded at 
t c . Meanwhile, according the proof above, we can easily see that 

m js (t)e- lx - x i-^ +m js+1 (t)e-\ x - x *> +1 W G J L 1 (M)nBV A (K) for all 1 < s < k, 

whose limit is 

m j .(t c )e- lx - x ^ t ^ + m Js+1 (t c )e^ x - x ^+ l( - t ^ 
-isgn(x-a: A (t c ))(i i .(t c )-i i . +1 (t c ))e-l*-^( t «)l 

as t — > t~, which leads to the conclusion. □ 

Corollary 5.5. The limit of a multipeakon u(-,t) for t —¥ t~ has a unique 
entropy weak extension which is a shockpeakon in the sense of H. Lundmark. 
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